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ABSTRACT 


By  means  of  an  automatic  digital  computer  solutions  were  ob¬ 
tained  for  a  number  of  second  order  servo  problems  with  backlash  in  the 
gear  train.  A  complete  range  of  solutions  were  obtained  with  variations 
in  load  friction,  motor  friction,  load  inertia,  motor  inertia,  damping  co¬ 
efficient  for  the  combined  operation ,  and  backlash  angle .  These  data 
are  presented  in  chart  form  suitable  for  the  servo  design  engineer  to 
estimate  performance  of  an  arbitrary  second  order  system  with  backlash. 

The  program  was  conducted  on  the  Control  Data  Corporation  Model 
1604  Computer,  and  this  pfogram  is  presented  with  explanations  of  the 
manner  in  which  the  differential  equations  were  solved.  From  data 
generated  by  the  computer  phase  plane  plots  are  drawn  as  a  further 
illustration  of  the  results  of  this  thesis ,  and  a  sample  problem  is  in¬ 
cluded  to  show  application  of  the  design  charts. 

The  writer  wishes  to  express  his  appreciation  to  Dr.  George  J. 
Thaler  of  the  Electrical  Engineering  Department  and  to  Dr.  E.  J.  Stewart 
and  Mr.  E.  N.  Ward  of  the  U.  S.  Naval  Postgraduate  Computer  Center 
for  their  assistance  and  encouragement  in  completing  this  work. 
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1. 


Introduction 


In  the  engineering  analysis  and  application  of  servomechanisms 
(hereafter  called  servos)  which  include  gear  coupling ,  the  problem  of 
backlash  must  be  considered  and  is  often  a  source  of  trouble.  In  theory 
it  is  desirable  to  know  the  system  performance  of  a  servo  with  a  given 
distribution  of  inertia  and  friction  and  to  accurately  determine  the  size 
of  the  limit  cycle  for  any  given  set  of  variables.  In  practice,  knowledge 
of  the  performance  of  a  system  may  be  desirable  prior  to  the  construction 
of  a  mock-up  or  simulation  system . 

Existing  methods  fonthe  analytical  solution  of  the  generalized 
backlash  problem  are  laborious  and  complicated  and  lack  of  the  straight¬ 
forward  approach  needed  in  the  preliminary  design  of  a  servo .  In  this 
thesis  design  charts  are  presented  for  the  second  order  servo  with  back¬ 
lash  to  accommodate  an  arbitrary  set  of  system  variables . 

Previous  work  in  this  problem  has  been  done  at  the  U.  S.  Naval 
Postgraduate  School  by  Lutkenhouse  in  reference  a  wherein  the  backlash 
problem  was  solved  by  construction  of  the  phase  trajectory  in  the  graphical 
phase  plane.  This  thesis  extends  the  work  of  Lutkenhouse.  It  is  recog¬ 
nized  that  the  construction  of  a  phase  trajectory  for  a  given  problem 
requires  many  hours  of  tedious  calculation  and  graphical  construction. 
Thus,  in  the  collection  of  a  large  volume  of  data  required  for  the  con¬ 
struction  of  design  charts  for  the  second  order  servo  with  backlash  an 
automatic  digital  computer  was  employed  . 
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The  work  of  this  thesis  first  required  derivation  of  the  system 
equations,  then  a  computer  program  was  developed  to  adequately  des¬ 
cribe  the  equations.  The  burden  of  the  work  subsequently  fell  in  pro¬ 
cessing  sufficient  data  to  construct  design  charts  to  satisfy  the  require¬ 
ments  of  the  problem. 
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2. 


Equations  for  the  Servo  Motion 


in  the  construction  of  a  set  of  design  charts  it  is  usually  desirable 
to  provide  sufficient  flexibility  so  that  the  charts  may  be  used  for  any 
given  servo.  Thus,  equations  are  presented  in  Appendix  A  for  the  de¬ 
termination  of  the  system  parameters  necessary  to  enter  the  design  charts 
of  this  thesis  for  three  types  of  second  order  servos:  the  armature  control 
motor,  the  field  control  motor,  and  the  two  phase  motor. 

All  of  these  cases  behave  in  a  similar  manner  when  coupled  to  a 
load  through  a  gearbox  which  has  backlash.  Initially  there  is  a  period 
when  the  backlash  is  taken  up,  and  the  system  performs  as  a  combined 
system.  The  differential  equation  describing  this  motion  for  an  undamped 
natural  frequency  of  one  is 

0  +  30  +0 -1 

where  B  is  equal  to  twice  the  damping  coefficient,  £ 

At  such  point  in  the  cycle  of  motion  wherein  the  load  velocity 
exceeds  that  of  the  motor,  backlash  starts,  and  the  motor  is  not  directly 
coupled  to  the  load.  There  are  two  differential  equations  describing 
this  motion.  The  equation  for  the  coasting  load  is 


where  D  = 


0i. 


JZ 


The  equation  for  the  motor  is 


where  E  =  J.„  and  F  =  F.  ,  . 

M  M 
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This  motion  in  the  backlash  region  continues  until  the  total  amount  of 
backlash  is  taken  up,  and  at  that  point  gear  contact  is  made.  The  action 
of  the  servo  upon  contact  is  not  a  simple  one.  This  motion  can  follow  one 
of  several  types  of  contact  where  the  basic  laws  governing  the  motion  are 
the  law  of  conservation  of  momentum,  and  the  law  of  conservation  of 
energy  which  must  both  be  considered.  In  this  thesis  it  is  assumed  that 
the  contact  is  a  plastic  one  with  no  gear  bounce;  therefore,  there  is  energy 
loss  in  the  gearbox  due  to  the  plastic  contact  and  the  law  of  conservation 
of  momentum  was  applied.  Following  the  gear  contact  phase  the  system 
is  again  combined,  and  the  cycle  of  events  begins  again. 

The  characteristics  of  the  second  order  servo  with  backlash  is 
independent  of  the  undamped  natural  frequency  as  shown  by  the  develop¬ 
ment  in  Appendix  A.  Therefore,  all  of  the  problems  worked  in  this  thesis 
and  used  in  the  construction  of  design  charts  are  based  on  an  undamped 
natural  frequency  of  one  radian  per  second,  and  the  step  input  is  equal 
to  one  radian.  Units  of  radians  and  seconds  are  used  throughout. 

The  amount  of  backlash  used  in  the  program  was  0.3  radian,  and 
it  was  a  purpose  of  the  thesis  to  determine  the  effect  of  other  magnitudes 
of  backlash. 
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3. 


Solution  of  the  Equations 


The  differential  equations  for  the  second  order  servo  with  backlash 
were  solved  on  the  Control  Data  Corporation  Model  1604  Computer.  This 
computer  is  described  in  reference  b,  but  briefly,  it  is  an  all-transistorized, 
stored  program,  general-purpose  digital  computer  having  a  large  storage 
capacity  and  exceedingly  fast  computation  and  transfer  speeds.  Input  to 
the  computer  was  by  means  of  magnetic  tape  and  punched  paper  tape,  and 
the  output  data  was  recorded  on  magnetic  tape.  The  data  was  then  printed 
on  the  International  Business  Machine  717  high  speed  printer. 

Figure  1  shows  a  schematic  of  the  problem.  In  programming  for  the 
solution  of  the  equations  it  was  necessary  to  successively  compute  a 
point  for  the  combined  motion,  test  this  point  to  determine  if  the  load 
would  drift  away  from  the  motor,  then  if  so,  continue  the  cycle  into  the 
backlash  region.  A  second  test  was  then  introduced  to  determine  the  end 
of  the  backlash  region,  and  if  the  backlash  was  fully  taken  up  the  contact 
part  of  the  program  was  initiated.  This  complete  program  was  accomplished 
by  the  program  of  Table  I . 

Several  instructions  in  the  program  are  worthy  of  mention.  First, 
it  is  noted  that  a  subroutine  was  used  in  the  actual  solution  of  the  dif¬ 
ferential  equations.  This  subroutine  based  on  the  Runge-Kutta  Gill  numer¬ 
ical  integration  is  contained  in  a  library  tape  in  the  U.  S.  Naval  Post¬ 
graduate  Computer  Center,  and  programming  instructions  are  contained 
in  reference  c.  The  program  was  assembled  in  accordance  with  the  ”AR" 
assembly  routine  of  references  d  and  e  which  also  contain  a  full  des¬ 
cription  of  each  instruction  used  in  this  program. 
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In  the  Runge-Kutta  Gill  numerical  integration  four  successive 
approximations  are  used  in  the  computation  oi  a  point  (in  this  thesis 
points  are  CL  01  second  apart) .  However,  it  was  not  desired  to  print  out 
each  point,  so  a  count  routine  was  inserted  in  the  program  to  print  each 
eighth  point.  Thus  in  the  actual  printout,  an  example  of  which  is  con- 
tained  in  Table  II,  test  points  are  at  intervals  of  .08  seconds.  The 
printout  format  shows  the  number  of  the  point  in  column  1 ,  the  time  in 
column  1,  the  motor  speed  in  column  3,  the  motor  position  in  column  4, 
the  load  speed  in  column  5,  and  the  load  position  in  column  6.  The 
number  following  the  figure  printed  represents  a  decimal  exponent  to 
locate  the  decimal  point. 

The  remaining  item  of  the  program  requiring  explanation  is  the 
portion  of  the  program  named  switch.  This  routine  is  inserted  to  accom¬ 
modate  the  various  changes  in  sign  caused  by  the  backlash  non-linearity, 

A  close  study  of  this  sub-routine  will  show  that  it  was  not  only  necessary 
to  change  the  sign  of  some  of  the  instructions,  but  it  was  also  necessary 
to  change  some  of  the  jump  instructions  caused  by  multiple  sign  changes 
in  the  backlash  test. 

With  the  basic  program  of  Table  I  stored  in  the  computer,  it  is 
merely  necessary  to  insert  the  desired  system  parameters  for  the  parti¬ 
cular  problem  to  be  solved  into  the  computer  locations  beginning  at  20156 
in  order  to  obtain  the  solution  for  any  desired  servo  problem  with  backlash. 
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4 .  Results 

A  total  of  120  runs  for  the  second  order  servo  with  backlash  were 
made  and  the  limit  cycle  lor  these  problems  were  determined.  Each  ran 
required  approximately  three  minutes  on  the  computer  and  yielded  one 
test  point  suitable  for  presenting  in  this  report.  These  test  points  are 
plotted  in  Figures  2  through  11,  Figures  12  through  15  are  phase  trajec¬ 
tories  for  a  number  of  selected  runs  plotted  from  an  actual  printout. 

Figure  2  is  a  plot  of  limit  cycle  versus  friction  ratio  for  a  damping 
coefficient  of  0.1  and  inertia  ratios  of  1/9,  1/4,  1,  4  and  9. 

Figure  3  is  a  plot  of  limit  cycle  versus  friction  ratio  for  the  load 
inertia  equal  to  the  motor  inertia  and  for  damping  coefficients  of  0.1, 
0.2,  0.3,  0.4,  0.5,  and  0,6. 

Figure  4  is  a  plot  of  limit  cycle  versus  backlash  angle  for  load 
inertia  equal  to  motor  inertia  and  for  a  friction  ratio  of  0.2.  This  plot 
shows  a  linear  variation  for  damping  coefficients  of  0.1,  0.2,  0.3,  0.4, 

0 . 5  ,  and  0.6. 

Figure  5  is  a  plot  of  limit  cycle  versus  backlash  angle  for  load 
inertia  equal  to  motor  inertia  and  for  a  damping  coefficient  equal  to  0.1, 
This  line  represents  all  friction  ratios  since  the  -  J  variation  of 
Figure  2  is  a  non-varying  straight  line.  For  a  damping  coefficient  other 
than  0,1  the  curve  would  spread  to  form  a  pencil  of  lines. 

Figure  6  is  a  plot  of  limit  cycle  versus  backlash  angle  for  a 
damping  coefficient  of  0.1  and  for  a  friction  ratio  of  0.2.  Variations 
of  the  curves  of  this  show  inertia  ratios  of  1/4,  1  and  4. 
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Figures  7  through  11  are  plots  of  limit  cycle  versus  inertia  ratio;: 
for  friction  ratios  of  0,1,  0.2,  0.4,  0.6,  and  1.0.  Each  plot  shows 
curves  for  variations  of  damping  coefficient  of  0.1  through  0.6. 

Figures  12  through  15  are  actual  plots  of  the  phase  trajectories  as 
plotted  from  the  print-out  sheets .  From  these  plots  can  be  determined  the 
output  speed  and  position  with  the  corresponding  time. 

Figure  12  shows  a  case  with  friction  ratio  of  0.2  at  an  inertia 
ratio  of  .25  and  a  damping  coefficient  of  0.1.  Figure  13  is  at  the  same 
friction  ratio  of  0 . 2  but  shows  an  inertia  ratio  of  4  and  a  damping  co¬ 
efficient  of  0.6.  Figure  14'and  15  illustrates  the  effect  of  increasing 
the  damping  coefficient  from  0.1  to  0.2  with  load  inertia  equal  to  motor 
inertia  and  at  a  friction  ratio  of  1.0. 
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Discussion 


5. 

After  the  accumulation  of  the  data  it  was  desirable  to  present  the 
information  in  the  most  useful  and  practical  form.  The  size  of  the  limit 
cycle  as  shown  in  the  figures  of  this  thesis  was  previously  defined  as 
deviation  from  the  equilibrium  or  average  value  rather  than  peak  to  peak 
variation.  The  data  were  plotted  in  various  forms  in  an  attempt  to  find 
the  most  useful  form.  Figure  2  is  one  such  presentation  where  the  data 
is  plotted  for  constant  damping  coefficient  versus  friction  ratio.  Figure 
2  was  not  found  to  be  the  optimum  form  because  in  subsequent  plots  at 
different  values  of  damping’coefficient  there  was  appreciable  crossing 
and  recrossing  of  lines  of  constant  inertia  ratio,  and  it  was  found  that 
interpolation  between  these  lines  was  difficult.  Further,  there  was  not 
a  wide  variation  of  limit  cycle  with  friction  ratio  as  evidenced  by  the 
level  curves  of  Figure  2.  In  fact,  for  the  particular  case  of  balanced 
load  and  motor  inertia  and  damping  coefficient  there  was  no  variation 
in  limit  cycle  with  friction  ratio.  This  plot  is  shown  as  a  logarithmic 
plot  because  it  was  found  that  the  digital  program  required  a  number  of 
modifications  to  satisfy  the  zero  friction  case.  A  single  case  was  run  with 
zero  friction  ratio ,  and  the  limit  cycle  was  found  to  be  substantially  equal 
to  that  for  a  friction  ratio  of  0.04  which  is  used  in  this  thesis  as  repre¬ 
sentative  of  the  zero  friction  ratio  case . 

Figure  3  is  a  composite  plot  showing  all  values  of  damping  co¬ 
efficient  up  to  0.6.  Figures  4,5,  and  6  are  presented  to  show  that  the 
size  of  the  limit  cycle  is  directly  proportional  to  the  amount  of  backlash . 


Several  runs  were  made  with  changes  in  only  the  backlash  to  conclusively 
prove  this  characteristic .  Figure  4  shows  the  variation  with  damping  co 
efficient  with  all  other  variables  constant.  Figures  5  and  6  are  similar 
plots  with  variations  in  other  factors ,  This  linear  variation  of  limit  cycle 
with  backlash  is  significant  since  it  means  that  design  charts  can  be 
drawn  for  any  constant  amount  of  backlash  and  variation  from  this  value 
is  in  direct  proportion.  In  this  thesis  all  charts  are  for  a  constant  amount 
of  backlash  equal  to  0.3  except  for  those  points  in  Figures  4,5,  and  6 
where  backlash  in  a  variable . 

Figures  7  through  11  hre  the  design  charts  found  to  be  the  most 
useful  form  for  design  estimations.  In  these  charts  the  size  of  the  limit 
cycle  is  plotted  versus  inertia  ratio.  The  logarithmic  scale  was  chosen 
for  two  reasons:  first,  it  avoided  the  zero  cases  found  at  the  extreme 
ends  of  the  scales  and  for  which  a  separate  program  would  have  been 
required,  and  secondly,  this  scale  gives  a  spread  to  the  curves  in  the 
region  of  inertia  ratio  equal  to  one  which  appears  desirable. 

The  design  charts  presented  herein  are  most  useful  for  the  range  of 
variables  encompassed  by  the  test  points.  While  it  is  believed  that  some 
estimates  can  be  made  from  these  charts  for  systems  with  very  large  and 
very  small  inertia  ratios  and  beyond  the  limits  of  the  curves  shown,  best 
results  are  obtained  for  values  where  interpolation  between  the  lines  of 
constant  damping  coefficient  is  feasible.  Appendix  B  gives  a  sample 
problem  using  the  design  charts .  Appendix  C  gives  instructions  in  the 
use  of  the  design  charts . 
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Conclusions 


6 . 

From  the  results  obtained  in  the  accomplishment  of  this  thesis 
it  is  concluded  that  the  size  of  the  limit  cycle  of  a  second  order  servo 
with  backlash  can  be  represented  by  design  charts ,  one  form  of  which 
are  presented  in  Figures  7  through  11  of  this  thesis. 

7 .  Recommendations 

Inasmuch  as  it  was  necessary  to  complete  this  thesis  during  a 
Summer  period,  time  was  of  essence,  and  the  scope  of  the  program  was 
definitely  limited.  Therefore,  the  following  recommendations  are  sub¬ 
mitted  as  a  reasonable  extention  of  the  program. 

1 .  Rewrite  the  digital  program  to  include  a  test  to  determine  the 
actual  size  of  limit  cycle  for  display  on  the  console  with  no  printout 
required.  This  will  permit  longer  runs  with  no  time  required  for  printout 
of  data.  Printout  capability  should  still  be  retained  by  means  of  a  Jump 
Key, 

2.  Extend  the  program  to  include  the  zero  friction  and  zero  inertia 
cases , 

3.  Extend  the  program  to  include  a  larger  range  of  inertia  ratios. 

4.  Reprogram  the  contact  phase  to  include  the  case  of  elastic 
contact  and  to  include  the  "lossy"  case  for  all  types  of  contact  between 
elastic  contact  and  plastic  contact. 
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9  „  Tables 

TABLE  I 

COMPUTER  PROGRAM 


REM  BACKLASH  PROBLEM  PLASTIC 

CONTACT  RUN  3 
ORG  20000 


20000 

75  4  20022 

START 

SLJ 

4 

INPUT  1 

REMARKS  IN  THIS 

00  0  00000 

0 

0 

0 

COLUMN 

20001 

75  4  20063 

COMBINED 

SLJ 

4 

RESET 

00  0  00000 

0 

0 

0 

20002 

75  4  60200 

SLJ 

4 

RUNGE 

SET  UP  GILL 

50  0  00000 

ENI 

0 

0 

(Procedure  to  solve 

20003 

00  0  20023 

0 

0 

TABLE  1 

differential 

00  0  20034 

0 

0 

DERC 

equations) 

20004 

75  4  60201 

POINTC 

SLJ 

4 

RUNGE  +1 

00  0  00000 

0 

0 

0 

20005 

75  4  60200 

LOAD  1 

SLJ 

4 

RUNGE 

SET  UP  GILL 

50  0  00000 

ENI 

0 

0 

2000G 

00  0  20040 

0 

0 

TABLE  L 

00  0  20051 

0 

0 

DERL 

20007 

75  4  60201 

POINTL 

SLJ 

4 

RUNGE  +1 

00  0  00000 

0 

0 

0 

20010 

75  4  20130 

SLJ 

4 

COUNT 

50  0  00000 

ENI 

0 

0 

20011 

75  0  20061 

SLJ 

0 

TEST  1 

50  0  00000 

ENI 

0 

0 

20012 

75  4  60200 

MOTOR 

SLJ 

4 

RUNGE 

SET  UP  GILL 

50  0  00000 

ENI 

0 

0 

20013 

00  0  20023 

0 

0 

TABLE  1 

00  0  20054 

0 

0 

DERM 

20014 

75  4  60201 

POINTM 

SLJ 

4 

RUNGE  +1 

00  0  00000 

0 

0 

0 

20015 

75  4  60200 

LOAD  2 

SLJ 

4 

RUNGE 

SET  UP  GILL 

50  0  00000 

ENI 

0 

0 

2001.G 

00  0  20040 

0 

0 

TABLE  L 

00  0  20051 

0 

0 

DERL 

20017 

75  4  60201 

POINT2L 

SLJ 

4 

RUNGE  +1 

00  0  00000 

0 

0 

0 

20020 

75  4  20130 

SLJ 

4 

COUNT 

50  0  00000 

ENI 

0 

0 

20021 

75  0  20073 

SLJ 

0 

TEST  2 

00  0  00000 

0 

0 

0 

13 


20022 

20023 

20024 

20025 

20026 

20027 

20030 

20031 

20032 

20033 

20034 

20035 

20036 

20037 

20040 

20041 

20042 

20043 

20044 

20045 

20046 

20047 
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75 

0 

00000 

INPUT  1 

SLJ 

0 

0 

Provision  for  changin  g 

75 

l 

20022 

SLJ 

0 

INPUT1 

Program  if  desired 

00 

0 

00000 

TABLE  1 

DEC 

2 

Table  for  gill  combined 

00 

0 

00002 

17 

7 

15075 

DEC 

.01 

34 

1 

21727 

00 

0 

00000 

T 

DEC 

0 

00 

0 

00000 

00 

0 

00000 

UDOT 

DEC 

0 

00 

0 

00000 

00 

0 

00000 

U 

DEC 

0 

00 

0 

00000 

00 

0 

00000 

QU 

DEC 

0 

00 

0 

00000 

00 

0 

00000 

THETADOT 

DEC 

0 

00 

0 

00000 

00 

0 

00000 

THETA 

DEC 

0 

00 

0 

00000 

00 

0 

00000 

QTHETA 

DEC 

0 

00 

0 

00000 

13 

0 

20156 

DERC 

LAC 

0 

B 

Derivative  for  combined 

32 

0 

20027 

FMU 

0 

U 

31 

0 

20032 

FSB 

0 

THETA 

30 

0 

20164 

FAD 

0 

ONE 

20 

0 

20026 

STA 

0 

UDOT 

12 

0 

20027 

LDA 

0 

U 

20 

0 

20031 

STA 

0 

THETADOT 

75 

0 

60202 

SLJ 

0 

RUNGE+2 

00 

0 

00000 

TABLEL 

DEC 

2 

Table  for  gill  load 

00 

0 

00002 

17 

7 

15075 

DEC 

.01 

34 

1 

21727 

00 

0 

00000 

TL 

DEC 

0 

00 

0 

00000 

00 

0 

00000 

WDOT 

DEC 

0 

00 

0 

00000 

00 

0 

00000 

W 

DEC 

0 

00 

0 

00000 

00 

0 

00000 

QW 

DEC 

0 

00 

0 

00000 

00 

0 

00000 

THETALD 

DEC 

0 

00 

0 

00000 

00 

0 

00000 

THETAL 

DEC 

0 

00 

0 

00000 

14 


20050 

20051 

20052 

20053 

20054 

20055 

20056 

20057 

20080 

20061 

20062 

20063 

20064 

20065 

20066 

20067 

20070 

20071 

20072 

20073 

20074 

20075 
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00  0  00000 
00  0  00000 

QTHETAL 

DEC 

0 

13  0  20157 

DERL 

LAC 

0 

D  Derivative  for  load 

32  0  20044 

FMU 

0 

W 

20  0  20043 

STA 

0 

WDOT 

12  0  20044 

LDA 

0 

W 

20  0  20046 

STA 

0 

IEETALB 

75  0  60202 

SLJ 

0 

RUNGE  +2 

13  0  20161 

DERM 

LAC 

0 

F  Derivative  for  motor 

32  0  20027 

FMU 

0 

u 

31  0  20047 

FSB 

0 

THETAL 

30  0  20164 

FAD 

0 

ONE 

33  0  20160 

FDV 

0 

E 

20  0  20026 

STA 

0 

UDOT 

12  0  20027 

LDA 

0 

u 

20  0  20031 

STA 

0 

THETADOT 

75  0  60202 

SLJ 

0 

RUNGE  +2 

50  0  00000 

EN1 

0 

0 

12  0  20046 

TEST  1 

LDA 

0 

T HEIALB  Test  to  determine 

31  0  20031 

FSB 

0 

THETADOT  Start  of  backlash 

22  2  20012 

JUMP1 

AJP 

2 

MOTOR 

75  0  20001 

SLJ 

0 

COMBINED 

75  0  00000 

RESET 

SLJ 

0 

0 

12  0  20025 

LDA 

0 

T  Changes  load  initial 

20  0  20042 

STA 

0 

TL  conditions  back  to 

12  0  20026 

LDA 

0 

UDOT  motor  initial  conditions 

20  0  20043 

STA 

0 

WDOT  if  separation  has  not 

12  0  20027 

LDA 

0 

U  occurred 

20  0  20044 

STA 

0 

W 

12  0  20031 

LDA 

0 

THETADOT 

20  0  20046 

STA 

0 

THETALD 

12  0  20032 

LDA 

0 

THETA 

20  0  20047 

STA 

0 

THETAL 

10  0  00000 

ENA 

0 

0 

20  0  20045 

STA 

0 

QW 

20  0  20050 

STA 

0 

QTHETAL 

75  0  20063 

SLJ 

0 

RESET 

50  0  00000 

ENI 

0 

0 

12  0  20047 

TEST2 

LDA 

0  ' 

THETAL  Test  to  determine  if 

31  0  20032 

FSB 

0 

THETA  backlash  has  been 

31  0  20166 

BKTEST2 

FSB 

0 

BACKLASH  taken  up 

50  0  00000 

ENI 

0 

0 

22  2  20137 

JUMP2 

AJP 

2 

CONTACT 

75  0  20012 

SLJ 

15 

0 

MOTOR 

20076 

20077 

20100 

20101 

20102 

20103 

20104 

20105 

20106 

20107 

20110 

20111 

20112 

20113 

20114 

20115 

20116 

20117 

20120 

20121 

20122 


TABLE  I  -  Continued 


36  0  20117 

SWITCH 

SSK 

75  0  20100 

SLJ 

75  0  20107 

SLJ 

50  0  00000 

ENI 

12  0  20062 

LEFT 

LDA 

14  0  20116 

ADD 

20  0  20062 

STA 

12  0  20075 

LDA 

14  0  20116 

ADD 

20  0  20075 

STA 

12  0  20074 

LDA 

15  0  20172 

SUB 

20  0  20074 

STA 

12  0  20117 

LDA 

05  0  00001 

ALS 

20  0  20117 

STA 

75  0  20001 

SLJ 

50  0  00000 

ENI 

12  0  20062 

RIGHT 

LDA 

15  0  20116 

SUB 

20  0  20062 

STA 

12  0  20075 

LDA 

15  0  20116 

SUB 

20  0  20075 

STA 

12  0  20074 

LDA 

14  0  20172 

ADD 

20  0  20074 

STA 

12  0  20117 

LDA 

05  0  00001 

ALS 

20  0  20117 

STA 

75  0  20001 

SLJ 

50  0  00000 

ENI 

00  1  00000 

FACTOR 

OCT 

00  0  00000 

25  2  52525 

ALTERNAT 

OCT 

25  2  52525 

75  0  00000 

PRINT 

SLJ 

12  0  20025 

LDA 

20  0  20252 

STA 

12  0  20027 

LDA 

20  0  20253 

STA 

12  0  20032 

LDA 

0  ALTERNAT  Plan  for  changing 

0  LEFT  sign  of  backlash 

0  RIGHT  each  cycle 

0  0 

0  JUMP  1 

0  FACTOR 

0  JUMP  1 

0  JUMP  2 

0  FACTOR 

0  JUMP  2 

0  BKTEST2 

0  FAC  TOR  2 

0  BKTEST2 

0  ALTERNAT 

0  1 

0  ALTERNAT 

0  COMBINED 

0  0 

0  JUMP1 

0  FACTOR 

0  JUMP1 

0  JUMP2 

0  FACTOR 

0  JUMP2 

0  BKTEST2 

0  FACTOR  2 

0  BKTEST2 

0  ALTERNAT 

0  1 

0  ALTERNAT 

0  COMBINED 

0  0 

0010000000000000 

2525252525252525 

0  0  Print  routine 

0  T 

0  BUF 

0  U  Motor  speed 

0  BUF+1 

0  THETA  Motor  position 
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20123 

20124 

20125 

20126 

20127 

20130 

20131 

20132 

20133 

20134 

20135 

20136 

20137 

20140 

20141 

20142 

20143 

20144 

20145 

20146 

20147 
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20  0  20254 

STA 

0 

12  0  20044 

LDA 

0 

20  0  20255 

STA 

0 

12  0  20047 

LDA 

0 

20  0  20256 

STA 

0 

75  4  71000 

SLJ 

4 

01  0  20252 

01 

0 

06  0  00001 

06 

0 

72  0  20126 

RAO 

0 

75  0  20120 

SLJ 

0 

75  0  00000 

COUNT 

SLJ 

0 

12  0  20170 

LDA 

0 

14  0  20171 

ADD 

0 

15  0  20173 

SUB 

0 

22  0  20134 

AJP 

0 

14  0  20173 

ADD 

0 

20  0  20170 

STA 

0 

75  0  20130 

SLJ 

0 

75  4  20120 

OK2PRINT 

SLJ 

4 

50  0  00000 

ENI 

0 

10  0  00000 

ENA 

0 

20  0  20170 

STA 

0 

75  0  20130 

SLJ 

0 

50  0  00000 

ENI 

0 

12  0  20162 

CONTACT 

LDA 

0 

32  0  20027 

FMU 

0 

20  0  20165 

STA 

0 

12  0  20163 

LDA 

0 

32  0  20044 

FMU 

0 

30  0  20165 

FAD 

0 

20  0  20031 

STA 

0 

20  0  20027 

STA 

0 

20  0  20044 

STA 

0 

20  0  20046 

STA 

0 

13  0  20156 

LAC 

0 

32  0  20027 

FMU 

0 

31  0  20032 

FSB 

0 

30  0  20164 

FAD 

0 

20  0  20026 

STA 

0 

13  0  20157 

LAC 

0 

32  0  20044 

FMU 

0 

20  0  20043 

STA 

0 

BUF+2 

W  Load  speed 

BUF+3 

THETAL  Load  position 
BUF+4 
DECO 
BUF 
1 

/-I 

PRINT 

0 

INDEX  Plan  for  printing 
INCRONE  only  every  tenth 
TEN  point  at  .08  sec 

OK2  PRINT  interval 
TEN 
INDEX 
COUNT 
PRINT 
0 
0 

INDEX 

COUNT 

0 

G  Apply  law  of 

U  momentum 

GU 

H 

W 

GU 

THETADOTPuts  resulting  values 
U  of  variables  in 

W  Table  1  and  L 

THETALD 
B 
U 

THETA 

ONE 

UDOT 

D 

W 

WDOT 
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20150 

20151 

20152 

20153 

20154 

20155 

20156 

20157 

20160 

20161 

20162 

20163 

20164 

20165 

20166 

20167 

20170 

20171 

20172 

20173 

20174 
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12  0  20047 

LDA 

0 

THETAL 

20  0  20032 

STA 

0 

THETA 

10  0  00000 

ENA 

0 

0 

20  0  20030 

STA 

0 

QU 

20  0  20033 

STA 

0 

QTHETA 

20  0  20045 

STA 

0 

QW 

20  0  20050 

STA 

0 

QTHETAL 

75  4  20120 

SLJ 

4 

PRINT 

50  0  00000 

ENI 

0 

0 

75  0  20076 

SLJ 

0 

SWITCH 

50  0  00000 

ENI 

0 

0 

50  0  00000 

17  7  66314 

B 

DEC 

.4 

B  equals  twice  the 

63  1  46314 

damping  coefficient 

17  7  54000 

D 

DEC 

.125 

D  equals  the  load 

00  0  00000 

* 

time  constant 

17  7  56314 

E 

DEC 

.2 

E  equals  motor  inertia 

63  1  46314 

17  7  64631 

F 

DEC 

.3 

F  equals  motor  friction 

46  3  14631 

17  7  56314 

G 

DEC 

.2 

G  equals  motor  inertia 

63  1  46314 

20  0  06314 

H 

DEC 

.8 

H  equals  load  inertia 

63  1  46314 

20  0  14000 

ONE 

DEC 

1. 

Integer 

00  0  00000 

00  0  00000 

GU 

DEC 

0 

Constant  initially  zero 

00  0  00000 

17  7  64631 

BACKLASH 

DEC 

.3 

Amount  of  backlash 

46  3  14631 

01  0  00000 

FACTOR2 

OCT 

0100000000000000 

00  0  00000 

00  0  00000 

INDEX 

OCT 

0 

00  0  00000 

00  0  00000 

INCRONE 

OCT 

1 

Used  in  count 

00  0  00001 

00  0  00000 

BLANK 

OCT 

0 

00  0  00000 

00  0  00000 

TEN 

OCT 

10 

Used  in  count 

00  0  00010 

17  7  15075 

1 

DEC 

.01 

Beginning  of  floating¬ 

34  1  21727 

point  conversion  table 
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20175 

20176 

20177 

20200 

20201 

20202 

20203 

20204 

20205 

20206 

20207 

20210 

20211 

20212 

20213 

20214 

20215 

20216 

20217 

20220 

20221 

20222 
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17  7 
34  1 

25075 

21727 

2 

DEC 

.02 

17  7 
12  1 

27534 

72702 

3 

DEC 

.03 

17  7 
34  1 

35075 

21727 

4 

DEC 

.04 

17  7 
63  1 

36314 

46314 

5 

DEC 

.05 

17  7 
12  1 

37534 

72702 

6 

DEC 

.06 

17  7 
60  5 

44365 

07534 

7 

DEC 

.07 

17  7 
34  1 

45075 

21727 

8 

DEC 

.08 

17  7 
07  5 

45605 

34121 

9 

DEC 

.09 

17  7 
63  1 

46314 

46314 

10 

DEC 

.10 

17  7 

36  5 

47024 

60507 

11 

DEC 

.11 

17  7 
12  1 

47534 

72702 

12 

DEC 

.12 

17  7 
72  7 

54121 

02436 

13 

DEC 

.13 

17  7 
60  5 

54365 

07534 

14 

DEC 

.14 

17  7 

46  3 

54631 

14631 

15 

DEC 

.15 

17  7 
34  1 

55075 

21727 

16 

DECq 

.16 

17  7 
21  7 

55341 

27024 

17 

DEC 

.17 

17  7 
07  5 

55605 

34121 

18 

DEC 

.18 

17  7 

75  3 

56050 

41217 

19 

DEC 

.19 

17  7 
63  1 

56314 

46314 

20 

DEC 

.20 

17  7 
36  5 

57024 

60507 

22 

DEC 

.22 

17  7 
12  1 

57534 

72702 

24 

DEC 

.24 

17  7 
00  0 

64000 

00000 

25 

DEC 

.25 
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20223 

20224 

20225 

20226 

20227 

20230 

20231 

20232 

20233 

20234 

20235 

20236 

20237 

20240 

20241 

20242 

20243 

20244 

20245 

20246 

20247 
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17 

7 

64121 

26 

DEC 

.26 

72 

7 

02436 

17 

7 

64365 

28 

DEC 

.28 

60 

5 

07534 

17 

7 

64631 

30 

DEC 

.30 

46 

3 

14631 

17 

7 

65075 

32 

DEC 

.32 

34 

1 

21727 

17 

7 

65463 

35 

DEC 

.35 

14 

6 

31463 

17 

7 

65605 

36 

DEC 

.36 

07 

5 

34121 

17 

7 

66314 

40 

DEC 

.40 

63 

1 

46314 

17 

7 

67024 

44 

DEC 

.44 

36 

5 

60507 

17 

7 

67146 

45 

DEC 

.45 

31 

4 

63146 

17 

7 

67534 

48 

DEC 

.48 

12 

1 

72702 

20 

0 

04000 

50 

DEC 

.50 

00 

0 

00000 

20 

0 

04314 

55 

DEC 

.55 

63 

1 

46314 

20 

0 

04631 

60 

DEC 

.60 

46 

3 

14631 

20 

0 

05146 

65 

DEC 

.65 

31 

4 

63146 

20 

0 

05463 

70 

DEC 

.70 

14 

6 

31463 

20 

0 

06314 

80 

DEC 

.80 

63 

1 

46314 

20 

0 

07146 

90 

DEC 

.90 

31 

4 

63146 

20 

0 

14000 

100 

DEC 

1.00 

00 

0 

00000 

20 

0 

14631 

120 

DEC 

1.20 

46 

3 

14631 

20 

0 

15463 

140 

DEC 

1.40 

14 

6 

31463 

20 

0 

16314 

160 

DEC 

1.60 

63 

1 

46314 

20 
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20250 

20  0  17146 

180 

DEC 

1.80 

31  4  63146 

20251 

20  0  24000 

200 

DEC 

2.00 

00  0  00000 

20252 

00  0  00000 

BUF 

BSS 

6 

00  0  00000 

DECO 

EQU 

71000 

RUNGE 

EQU 

60200 

END 
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1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


TABLE  II 


TYPICAL  PRINTOUT  OF  DATA 


MOTOR 

SPEED 

77040630-  1 
15278417 
22518046 
29387916 
35856424 
41895488 
47480616 
52590944 
57209252 
61321964 
64919110 
67994283  ’ 
70544563 
72570424 
74075629 
75067101 
75554775 
75551450 
75072608 
73561932 
70007602 
64634564 
57670572 
49317406 
39753812 
29138106 
17610490 
52951014-  1 
-76981710-  1 
-21272070 
48719261 
44993108 
39850916 
34619795 
29335478 
24032905 
18746033 
13507650 


MOTOR 

POSITION 

30333708-  2 
12247439-  1 
27389548-  1 
48177680-  1 
74303153-  1 
10543340 
14121485 
18127574 
22522912 
27267569 
32320675 
37640704 
43185760 
48913847 
54783141 
60752245 
66780440 
72827914 
78855989 
84814975 
90570699 
95967717 
10086982  1 

10515798  1 

10872835  1 

11149055  1 

11336614  1 

11428723  1 

11419529  1 

11304005  1 

14276023  1 

14557350  1 

14896797  1 

15194727  1 

15450572  1 

15664046  1 

15835140  1 

15964112  1 


LOAD 

SPEED 

68901844-  1 

14489345 

21757968 

28660703 

35165619 

41244294 

46871885 

52027169 

56692568 

60854140 

64501553 

67628038 

70230320 

72308523 

73866071 

74909556 

75448598 

75495690 

75066030 

74326600 

73587038 

72854835 

72129917 

71412213 

70701649 

69998156 

69301663 

68612100 

67929398 

67253489 

48719261 

45474833 

40349320 

35131214 

29856282 

24559521 

19274963 

14035491 


LOAD 

POSITION 

30007982-  2 
12215856-  1 
27359124-  1 
48148568-  1 
74275496-  1 
10540733 
14119047 
18125317 
22520843 
27265695 
32319002 
37639237 
43184500 
48912797 
54782300 
60751613 
66780013 
72827688 
78855961 
84831950 
90748447 
96606073 
10240541  1 

10814705  1 

11383155  1 

11945950  1 

12503145  1 

13054795  1 

13600957  1 

14141684  1 

14276023  1 

14557542  1 

14896996  1 

15194932  1 

15450780  1 

15664257  1 

15835352  1 

15964323  1 


22 


flgtv  Eq£?  B^o<la5H 


/^/(J7C/£’ £- / 


33. 


APPENDIX  A 


DEVELOPMENT  OF  CHART  PARAMETERS 
The  purpose  of  this  appendix  is  to  present  a  procedure  whereby  the 
factors  required  to  enter  the  design  charts  of  this  thesis  may  be  computed 
for  a  given  servo  system.  Equations  for  the  field  control  motor,  the  arma¬ 
ture  control  motor,  and  the  two  phase  motor  are  developed.  The  application 
of  a  given  system  represented  by  a  differential  equation  of  the  form 


requires  a  transformation  to  the  form 

since  the  design  charts  are  based  on  an  undamped  natural  frequency  of  one. 
This  can  be  accomplished  by  a  time  scaling  of  the  problem.  Thus,  in  the 
first  equation  above  if  the  independent  variables  are  related  by  ' 

where  ,  tfg  dll  —  J!  , 

Ji  ^ 

and  $  =  d 

When  these  substitutions  are  made  there  results  the  desired  form  suitable 
for  use  with  the  design  charts.  From  this  ,  it  is  evident  that  the  damping 
coefficient,  ,  remains  unaltered  in  the  transformation,  and  the  design 
charts  ,  while  based  on  an  undamped  natural  frequency  of  one,  are  applicable 
to  other  frequencies  as  well. 

The  second  order  servo  can  be  represented  by  a  block  diagram  as 
follows . 
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Three  types  of  motors  are  considered,  the  armature  control  motor,  the 
field  control  motor,  and  the  two  phase  motor.  For  the  period  of  time  when 
the  gears  are  in  contact  the  system  is  a  simple  second  order  system.  The 
differential  equations  for  this  combined  motion  can  be  derived  by  equating 
the  load  torque  required  to  the  motor  torque  developed.  For  the  armature 
control  motor  this  results  in  the  following  equation. 


For  the  field  control  motor  if  the  field  time  constant  is  negligible  the 


equation  is  similar. 

In  like  manner  the  differential  equation  for  a  two  phase  motor  is 


ZZ  y~ZZ 

where 


Ki  =  motor  torque  constant 

K  ,  =  motor  generator  constant 

=  error  measurement  constant 
R  =  armature  resistance 
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f*  =  gear  constant 

T  ,  =  inertia  of  motor  referred  to  load  shaft 

M 

J  =  inertia  of  motor 

Jm 

J  =  inertia  of  load 

L 

fm  =  friction  constant  of  motor 

=  friction  constant  of  motor  referred  to  load  shaft 
F  =  friction  constant  of  load 

Fj  =  friction  of  load  plus  motor  friction  referred  to  load  shaft 
v  _  £3 


(from  characteristic  curves  of  two  phase  motor) 

(from  characteristic  curves  of  two  phase  motor) 


r/? 


K 


N 


9T 

Bn 


K, 


ar 
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Procedure  for  Armature  Control  Motor 


1. 

2. 

3. 

4. 

5. 

6. 

7. 

8. 

9. 


Find  K  from 


/<r= 


X/  X-7 


Find  J  referred  to  load  shaft . 
M 


Find 


y/7 


*4 


Find  Ft  from  /T_  =-  ^  y. 


Find  F^  from 
Find  4 
Find  JL/JM 


Find  FL/FTotal  = 


Use  6,7,  and  8  to  select  proper  design  chart  and  find  size  of  limit 
cycle  for  0.3  rad  backlash. 


10.  Find  size  of  limit  cycle  for  desired  backlash  from 


zC/aV/7"  A /a//,  C/ArAfJ  [  y 


Procedure  for  Field  Control  Motor 


1.  through  4.  Same  as  Armature  Control  Motor 

5 .  Find  4  ^ 

6.  Continue  as  Armature  Control  Motor 


Procedure  for  Two  Phase  Motor 


1.  Find  K  from  characteristic  curves  of  motor 
e 


2 .  Find  T ,  referred  to  load  chart 

M 

3 .  Find  Ojn  = 


4.  Find  F 


total 


5.  Find  K  from  characteristic  curves  of  motor, 
n 
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7.  Continue  as  Armature  Control  Motor. 
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APPENDIX  B 


SAMPLE  PROBLEM 

A  Elinco  FD-162  115  volt  D.  C.  Shunt  Motor,  1/125  Horsepower 

4,000  Rpm  is  coupled  through  a  gearbox  to  a  load.  The  constants  of  the 

proportional  error  servo  are  as  follows: 

Armature  resistance,  R  —  282  ohms 

Motor  Generator  constant,  K^O.233  volts/rad/sec 

Torque  Constant,  K«28.2  ounce  -  inch/ampere 

2 

Motor  Inertia,  J  ==  0.61  ounce -inch 
m 

-s 

Motor  Friction,  f  —.00064  ounce-inch/rad/sec 
m 

Error  measurement  constant,  Kg=10  volts/rad 
Gearbox  factor,  -^>=2 
Backlash  0 . 1  rad 

-6  2 

Load  Inertia  J  =16.4  x  10  slug  ft 

L 

Load  Friction  Fj^.002  ounce-inch/rad/sec 


Solution:  Using  the  procedure  for  the  armature  control  motor  in  appendix  A. 

1  //=  —  —  C^S.B)(/o) -  /r-r- 


2.  =  -4?r  = 


£>  / 


—  3. 2  x/o  C  Jx-c/G 


=  l  /  - =  1/  — -  &  £&^/sec 
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6. 


>  /-^b ±JlL  V-L-  j  1/^  __ 

>  5^vp/  6&.6*,oyc<zX'4*J  -/9 


/_ 

r%/,} 


x/O 


7. 


J2"  __  /£.4  x/& 

S.<^y  /&'<=> 


=  a 


8. 


^2 

/ZT577^ 


=•  .077 


9 .  From  Figure  7  representing  a  friction  ratio  of  .  1 ,  and  for  an 
inertia  ration  of  2.0  and  a  damping  coefficient  of  .19  the 
size  of  the  limit  cycle  for  0.3  radian  backlash  is  .25  radian. 

10.  The  size  of  the  limit  cycle  for  0.1  radian  backlash  is 
limit  cycle  =  -£rf-  ~  •  *#3  . 

0  3/  -  ■ 
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APPENDIX  C 


USE  OF  DESIGN  CHARTS 

Figures  7  through  11  are  presented  as  representing  design  charts 
for  the  second  order  servo  with  backlash.  These  charts  can  be  used  to 
rapidly  determine  the  size  of  the  limit  cycle  by  the  following  procedure . 

1.  If  the  system  variables  are  not  known,  compute  these  by 

the  methods  of  Appendix  A. 

2.  Choose  the  chart  representing  the  Load  Friction/Total  Friction 

desired. 

3.  Read  the  Size  of  the  Limit  Cycle  resulting  from  the  intersection 

of  the  abscissa  representing  inertia  ratio  and  the  desired 
damping  coefficient,  interpolating  if  necessary. 

4.  The  value  determined  represents  that  Limit  Cycle  resulting 

from  a  backlash  of  0.3  radians.  For  a  given  backlash  the 
limit  cycle  is  directly  proportional,  thus  divide  the  Limit 
Cycle  found  in  3  above  by  0.3  and  multiply  by  the  desired 
backlash. 


45 


